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ABSTRACT 

Solar wind fluctuations, such as magnetic field or velocity, show power law 
power spectra suggestive both of an inertial range of intermittent turbulence 
(with ~ —5/3 exponent) and at lower frequencies, of fluctuations of coronal ori- 
gin (with ~ —1 exponent). The ULYSSES spacecraft spent many months in 
the quiet fast solar wind above the Sun's polar coronal holes in a highly ordered 
magnetic field. We use statistical analysis methods such as the generalized struc- 
ture function (GSF) and extended self-similarity (ESS) to quantify the scaling 
of the moments of the probability density function of fluctuations in the mag- 
netic field. The GSFs give power law scaling in the range of the form 
(\y(t + r) — y{t)\ m ) ~ T^ m \ but ESS is required to reveal scaling in the inertial 
range, which is of the form (\y(t + r) — y(t)\ m ) ~ [g(T)]^ m \ We find that g(r) 
is independent of spacecraft position and g(r) ~ 7 --'°sio(At) _ rp^g "J- 1 " scaling 
fluctuates with radial spacecraft position. This confirms that, whereas the 
fluctuations are directly influenced by the corona, the inertial range fluctuations 
are consistent with locally evolving turbulence, but with an "envelope" g(r), 
which captures the formation of the quiet fast solar wind. 

Subject headings: magnetic fields - solar wind - turbulence 



1. Introduction 



The solar wind provides a uniqu e opportunity to perfo rm in-situ long interval studies 
of a high magnetic Reynolds number (IMatthaeus et al.ll2005l ). magnetohydrodynamic stellar 



1 Centre for Fusion, Space and Astrophysics, Department of Physics, University of Warwick, Coventry, 
CV4 7AL, UK 

2 UKAEA Culham Division, Culham Science Centre, Abingdon, Oxfordshire, 0X14 3DB, UK 



-2- 



wind. Magnetic field fluctuations i n the solar wind typically exhibit an ex tended power law 



region in the power spectrum (e.g, IColemanlll968l ; iBruno fc Carbond 120051 ) . from timescales 



of a few hours down to that characteristic of the ion dynamics. Solar wind acceleration and 
heating to form t h e fast solar wind is known to occur on open flux lines, i.e. over coronal 
holes (see lOfmanl (120051 ) for a recent review). The solar polar coronal holes thus provide 



an uninterrupted, spatially extended region in which we can study the accelerated fast solar 
wind. 

A key objective of the present work is to relate aspects of the ULYSSES solar wind 
measu rements described in the next section: their spectral power density; their intermit- 
tency (IBruno et al.ll2003l ); and their spatial location. Generalizing s omewhat, the solar wind 
spectral power density is observed to scale approximately as (IRuzmaikin et al.lll9 95b 



Goldstein et al 



1995ah_ at l ower frequencies (< 1 mHz); and as "/ 5//3 " (IRuzmaikin et al. 



19931 ; lliorbury et al. 1995a ). reminiscent of the inertial range of Kolmogorov (jKolmogorov 
19411 ). at higher frequencies (~ 10 mHz- 100 mHz). The frequency at which the transi- 



tion occurs (~ 1 mHz-10 mHz) between these two power laws is ob served to decline with 



increasing distance f rom the sun in the plane of the ecliptic (see iFeynman et al.l 11995 
Horbury et al.l Il996al ). This extension of the "J -5 / 3 " range to lower freq uency at greater 



distances can be interpreted as evidence for an active turbulent cascade (IGoldstein et al. 



1995al ; iHorbury &: Baloghl 119971 ) that evolves a growing inertial range as time passes in the 
outward propagating plasma. The compone nt is taken to reflect embedd ed coronal 

turbulence, convected with the solar wind plasma (IMatthaeus fc Goldstein! Il986l ). It is not 
yet certain whether, in addition, this coronal turbulence acts as the low frequency large- 
scale driver of the inertial range turbulence. Both the and "J -5 / 3 " fluctuations are 
often predominantly shear Alfvenic in character, that is incompressible and displaying cor- 
relation or anticor relation between perturbations of the magnetic field and of fluid velocity 
jSmith et al.lll995l) . 



The first observations initiated a debate as to whether the fluctuations are of solar coro- 
nal origin, simply pas sively advecting with the flow, or whether they are dominated by locally 
evolving turbulence (IGoldsteinl 120011 ) . The inertial range fluctuation s, and the crossover to 
"j-i" behavior, show secular variation with heliographic distance (IHorbury et al.l Il996al ) 
consistent with evolving, rather th an fully evolved, turbule nce. Helios data, in conjunction 



with ULYSSES, has been used by IGoldstein et al.l (jl995bl ). to show that the "/ 1 " region 



contains scales, which are too large to be produced in situ and are progressively eaten away 
by the smaller scale turbulence, which must therefore be active. The large scale magnetic 
structure of the corona also varies with both heliospheric latitude and solar cycle, and this is 
clearly manifested in the coherent structures and variation of wind speed that are observed 
(jPhillips et al.lll995l ). Power spectra are not sufficient to quantify fully the scaling properties 
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of fluctuations (e. g. IChapman et al.ll2005l ) 



Recent structure function analysis using WIND and ACE data, taken in the eclip- 
tic plane at 1AU, shows evid ence of scaling within the inertial range that is solar cycle 
dependent ( jKiyani et al.l 120071 ) . To unravel the interplay between the large scale coronal 
driver and the evolving inertial range turbulence, we therefore make use of ULYSSES polar 
passes. ULYSSES' out-of-ecliptic orbit enables the study of possible lati tudinal and ra- 
dial dependences of the solar wind (ISmith et al.lll995l ; iHorbury et al.lll996al ). In particular, 
ULYSSES spent many months in the quiet fast solar wind above the polar coronal holes 
at periods of both minimum and maximum solar activity. The passes during periods of 
minimum solar activity have been extensively studied because they allow observations of 
magnetohydrodynamic (MHD) fluctuations free from perturbations from large scale events 
such as coronal mass ejections (CMEs) or large scale stream structures. The instruments 
aboard ULYSSES relevant to studies of solar wind fluctuations are the Vector Helium Mag- 
netometer (Jet Pro pulsion Laboratory) and the Fluxgate Magnetometer (Imperial College) 
(Bal ogh et al. 1992). 3 lasma flow measurements are obtained from the SWOOPS experiment 
(jBame et al.lll992l ). and the time resolutions available range from seconds for the magnetic 
field measurements to 4 minute averages for the electron and ion velocities. Thus the mag- 
netometer data offers a sufficiently extensive dynamical range to explore t he inertial range. 
ULYS S ES has performed five polar passes and p r eviou s studies include ISmith fc Balogh 
(119951 ); IHorbury et al.l (Il996bl ); IHorbury fc Baloghl (120011 ) and references therein. 



In this paper, we seek to characterize the intermittency therefore we remove the mini- 
mum of outliers consistent with obtaining goo d representation of the tails of the probability 
density functions (PDFs) (jKiyani et al.ll2006l ). This is distinct from, but complimentary to, 
approaches that seek to eliminate the intermittency by removing significant fractions of the 
tails of the PDFs, in order to probe the remnant scaling (S. Haba l, private co mmunication, 
2007). We perform generalized structure function analysis (GSF) (jFrischlll995l ). on intervals 
of quiet solar wind as seen by ULYSSES. We show that wh ile GSF is not sufficient to reveal 



scalin g in the inertial range, extended self-similarity (ESS) (IBenzi et al.lll993l ; ICarbone et al. 



19961 ). successfully recovers self-similar behavior. Furthermore we examine the possible lat- 
itudinal and radial dependences of both the inertial and ranges, and conclude that 
the inertial range, unlike the regime, shows very little variation with spacecraft po- 
sition over the 60 days considered. We focus here on magnetic field fluctuations for our 
analysis rather than Elsasser variables, s i nce t hese involve velocity measurements at lower 
cadence; however see ISorriso-Valvo et al.l (120071 ). where local (in time) scaling properties are 
also considered. 
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2. Quantifying Scaling 



We first introduce the analysis tools that enable us to quantify s caling in a timeseries. 
Generally speaking, a time series ?/(£) exhibits scaling (jSornettell2004l ). if 



Here the angular brackets denote an ensemble average over t, implying an assumption of 
approximate stationarity In practice, we examine the data over a sufficiently large range of 
time intervals r in order to establish the power law dependence (see equation (0Q)); that is, 
the scaling exponent ({m). 

The poin t of contact between equation (p]) and fluid turbulence is the Taylor hypothesis 
(jTaylorlll938l ). where in the high speed flow the time interval r plays the role of a longitudinal 
lengthscale. However equation (pQ) expresses a generic scaling property, whi ch is also found , 
for example, in random fractals such as Brownian walk s and Lew flights (|Sornettd 120041 ) . 
Since the coronal magnetic carpet is known to be fractal (jSchrijver et al.lll992[) and the larg e 
scale coronal dynamics as seen in solar flare statistics exhibit scaling (jCrosby et al.lll993l ). 
one might anticipate that a propagating signature of coronal origin might also show scaling 
which could be captured by equation ((Tj). In practic e, we test equation (HD by computing the 



associated generalized structure functions or GSF (IBurlaga fc Kleinlll986l ; iRuzmaikin et al. 



19931 . Il995al ; iMarsch fc Tul Il996l ; iHorbury fc Baloghl 119971 . and references therein) : 



S m (T) = {\y(t + r)-y(t)n = {\5yn 



(2) 



For a perfectly self-similar process, such as the Kolmogorov cascade or a random fractal, 
C(m) depends linearly on m. Turbulent flows are however typically intermittent, that is, the 
dissipation does not occur unifor mly in space or time. The corresponding exponents ((m) 
are quadratic in m (IFrischl Il995l ). Nevertheless, S m scaling with r would be expected in 
uniform, fully developed turbulence in an infinite domain. 



In practice, both in numerical simula t ions (IMerrifield et al.ll2007l ) and in the laboratory 
(IDudson et al.l 120051 ; iDendy fc Chapman! 120061 ) the scaling in equation ([TJ is not always 
found. This may reflect the finite spatial domain, or that the turbulence is not fully evolved. 
We will also see that this is the case for the inertial range in quiet intervals of ULYSSES 
observations. However a weaker form of scaling, known as extended self-similarity (ESS) 
turns out to be applicable to our datasets in cases where equation is not. ESS proceeds 
by replacing r in equation ([T]) by an initially unknown generalized timescale g(r), such that 
formally 



S m (r) ~ [g(r)f m ~> 



(3) 
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It follows from equation ([3]) that 



(4) 



see for example iGrossmann et al.lll997l ; iPagel fc Baloghll200ll ). 



The measured vector magnetic field time series, B(i), is differenced for time lags r in 
the range 1 minute to 50 minutes, yielding a series 5yi(t, r) for its three components 



6y i (t,T) = B i (t + T)-B i {t) 



(5) 



where % denotes the component of B under consideration. Assuming time-stationarity, the t 
dependence in Syi(t,r) can be dropped: 



T 



(M m ) 



\6y l \ m P(Sy l ,T)d(6y l ) 



(6) 



where P{5y^ r) is the probability density function of Syi, < ■ > again denotes temporal 
averaging and m is a positive integer. The effect of outliers becomes increasingly apparent in 
the higher order structure functions, subject to statistical constraints, which typically limits 
consideration to m < 6. One of the main problem s associated with the use of structure 
functions is the nature of the limits in equation (| 6|) (jChapman et al.ll2005l ). For this reason 
a clipping technique developed by Kiyani et al. ( 20061 ) is applied, when indicated below, 
to condition the data. This involves removing an increasing fraction of the maximum and 
minimum y^ values and observing the effect on the scaling. Both raw datasets and conditioned 
datasets will be shown. 



3. The datasets 

In 1995, the ULYSSES spacecraft spent three months above the North polar coronal 
hole, in quiet fast solar wind. As this was close to a period of minimum solar activity, the 
magnetic topology of the Sun was relatively simple, free from transient events such as solar 
flares. The Sun's surface magnetic field was actually dipolar during this time, positive or 
outwards in the Nor thern hemisphere and negative or inwards in the Southern hemisphere 



f lForsvth et al.lll996h . 



Throughout this paper, we present results separately for each 10 day contiguous interval 
from day 180 to 239 of 1995, while ULYSSES was above the Northern coronal hole. Each 
dataset comprises approximately 13, 000 datapoints of one minute averaged measurements. 
This enables us to explore both the scaling properties of the inertial and ranges and 
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to test for radial and latitudinal dependencies. Over these 60 days, ULYSSES moved from 
1.7926 AU to 2.2043 AU heliospheric distance and through 73.76° to 77.03° heliographic 
latitude, with a peak at 80.22° on days 212 and 213. The successive time intervals are 
compared in order to identify any radial or latitudinal trend. The study is restricted to the 
radial (R), tangential (T) and normal (JV) magnetic field components and uses one minute 
averaged measurements in order to remove any possible sub-spacecraft spin artefacts since 
ULYSSES has a spin period of 12s. The RTN coordinate system corresponds to solar- 
ordered coordinates where R is the sun-ULYSSES axis, T is the cross product of R with 
the solar rotation axis, and N is the cross product of R and T, completing the right-handed 
system. 



4. ULYSSES observations and scaling 
4.1. Power Spectra 



The magnetic field power spe ctra are computed using the multitapering spectral analysis 
method (jPercival fc Waldenlll993l ). In Figure[T]we see that the power spectra show an inertial 
range with a Kolmogorov-like behavior at higher frequencies and a characteristic flattening of 
the spectra at lower frequen cies The existence of this regime is well-known in many physical 



proce sses ( IBak et al.l Il987l ). including the interplanetary magnetic field (IMatthaeus et al. 



20071 ). The power spectra reveal power law scaling, but give no information on intermittency 
or on whether the turbulent cascade is active; for this we turn to the associated GSFs. 



4.2. Generalized Structure Functions 

We first summarise the scaling behavior seen in the different time intervals in Figure [21 
where we pl ot versus r for th e data conditioned by clipping 0.1% outliers following the 



technique of iKiyani et al.l (120061 ) . The existence of two distinct scaling regions is clear. For 
the small r region we see the inertial range with scaling exponents Cs( m )> whereas the large 
t region corresponds to the range with scaling exponents (^(m). Figure [3] shows that 

linear regression applied to the third order structure functions for both raw and conditioned 
data yields power-law scaling in the range. However this does not appear to be the 

case for the inertial range, where the S m clearly curve for r < 30 minutes. The ( exponent 
numerical values for the raw and conditioned GSF data differ by no more than ~ 6% on 
average, with the higher order moments showing stronger variations, as expected. 



Figures [2] and [3] suggest that a single function g(r) may be common to all the time 
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intervals considered. We test this conjecture in Figure HI by normalizing all the S3 log plots 
to their value at r = 30min., close to the centre of the r range. The curves overlay quite 
closely, and there is no secular latitudinal or radial dependence of g{r). The displacements 
of the curves in Figures [2] and [3] arise from the data: as ULYSSES moves, the GSF and ESS 
plots shift in a relatively ordered way. In Figure [5] we examine this further by showing the 
variations of a single point on each line, S${t = 30) as a function of time and therefore of 
increasing radial distance, for all B field components. S3 (r = 30) for all three components 
exhibits a secular trend and decreases with time as heliographic range increases. There is no 
significant latitudinal variation, although there is a flattening of the slope after ~ day 210 
(or interval 4), when ULYSSES passes over the solar pole and the heliographic latitude starts 
to decrease again. This is more apparent in the iV and T comp onents, the R componen t 
seeming relatively unaffected. This is consistent with the work of iGoldstein et al.l (Il995bl ). 
who also observed a stronger radial, rather than latitudinal, dependence of the turbulent 
properties. 



4.3. Extended self-similarity 

Let us now test more precisely for S m ~ [<7(7")] within the inertial range by applying 
ESS to the data. We first apply the technique to the entire r range of 2 — 49 min., with 
results shown in Figures El and [3 

Power law scaling is recovered for the entire r range considered. ESS also seems to 
extend the inertial range scaling region and, apart from a few exceptional cases, it is generally 
difficult to distinguish clearly between the inertial range and ranges from these plots. 

Figure shows a reasonable linear fit to the whole r range but we see a small but systematic 
displacement of the data from the fitted line at higher frequencies. Figure [7] shows the ESS 
exponents £(2)/£(3), obtained from the gradients of the linear fits in Figure El which exhibit 
no dependence on the time interval considered. 

However the GSF analysis has given prior indication of the r at which the transition from 
one regime to the other occurs. It is therefore of interest to apply ESS analysis separately to 
the two regions - inertial range and - identified above. As we discuss below, Figures 

IH1 to [TOl then demonstrate that the range is the dominant source of variation, in clear 

contrast to the inertial range. Figure [H] repeats the ESS analysis for the inertial range with 
r = 2 — 14 min., and Figure [9] for the "J 1 " range with r = 26 — 49 min. 

The difference between the inertial range and fits can be seen in Figure [101 where 

the gradients of the fits are shown. Importantly, Figure [TD] demonstrates that the dominant 
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contribution to the scatter in the corresponding global plot, Figure [5], arises from the "/ x " 
region at larger r. This aligns with what was found with the GSF analysis, see Figured 



The GSF plots in Figure [2] can be fitted using a quadratic fit of the form 
log w S m (r) ~ a(m)(log 10 T) 2 + (3(m)(log 10 T) + 7 



(7) 



where 7 contains the secular variation in S3 seen in Figure Equation (j7j) can be rewritten 
in the form 



T 



-a(m)(logio(r)+X) 



-log 10 (\T)]a(m) 



where A = 10 A = io~ 7 - 575±0 - 246 is found by taking an average over fits for all the different 
time intervals. The minus sign is necessary if we insist on positive scaling of exponents a(m) 
as a is negative and (3 is positive. In Figure [111 we show A for all components and time 
intervals. The recovery of a power law behavior in the inertial range can be seen in Figure 
[T2lby plotting the GSF versus the g{r) defined in equation ([8]) rather than r, where the form 
of g(r) above is normalized such that Cs(3) = 1- 

Finally let us again consider the behavior of the (s exponents for m = 1 — 6 for all 
components and time intervals. These can be fitted to the multifractal model for fully 
developed turbulence proposed by iMeneveau fc Sreenivasanl (119871 ) (p- model). A simple 
example of this is to consider an eddy decaying into two smaller ones. The parent transfers a 
fraction p of its energy to one of them and a fraction (1 —p) to the other one (by convention 
p > 1/2). The values of ((m) are given by 



C(m) = 1 - log 2 (p m / 3 + (1 - p) m ") 



(9) 



which reduces to Kolmogorov scaling, i.e ((m) = m/3, for p = 1/2. For the data considered 
here, p-model fits of the (s exponents for the different components during the first time 
interval give p R = 0.79 ± 0.01, pr = 0.85 ± 0.01 and p/y = 0.84 ± 0.01. These values are 
consistent with previo us work by iHorbury fc Baloghl (jl997l ). who found p ~ 0.8, and by 



Paeel fc Baloghl faOOA 



5. Conclusions 

In this paper we have analysed ULYSSES quiet fast polar solar wind magnetic field 
measurements to study the evolving turbulence. We quantify the scaling behavior of both 
the inertial range (i.e. power spectrum ~ /~ 5/ ' 3 ) and the lower frequency range present 

in the solar wind. Six contiguous intervals of ten days, over which approximate stationarity 
can be assumed, were studied using the Taylor hypothesis to relate temporal and spatial 
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scales. The scaling that we establish is "macroscopic" in the sense that it is obtained over 
these time intervals of ten days. One can also consider local (in time) scaling properties 



( ISorriso-Valvo et al.l 120071 ). We use generalized structure functions (GSF) and extended 



self-similarity (ESS) to quantify statistical scaling, and find that: 

• GSF is sufficient to reveal power law scaling in the low frequency range of the 
form {\y(t + r) — y(t)\ m ) ~ r^ m \ but ESS is necessary to reveal scaling in the inertial 
range. This implies a scaling of the form (\y(t + r) — y(t)\ m ) ~ [^(t)] over the 
inertial range. 

• The range scaling varies in a non secular way with spacecraft position as found 



previously (IHorbury et al.l Il995bl ) . This is consistent with a coronal origin for the 
"J- 1 " scaling. 

In the inertial range, comparisons of the third order structure function S3 for the 
different time intervals show that g{r) is independent of spacecraft position, although 
an ordered trend in the absolute value of S3 with increasing radial distance is observed. 



A good fit to the inertial range is g(r) ~ r i °9io( Ar ) ) where A = 10" 



7.575±0.246 



• The exponents found for the inertial range, normalized such that Cs(3) = 1, are fitted 
by a p-model with p R = 0.79 ± 0.01, p T = 0.85 ± 0.01 and p N = 0.84 ± 0.01. This 
implies a higher degree of intermittency in the normal components of the magnetic 
field than in the radial component. 

Our results clearly differentiate between the dynamics of the fluctuations seen in the 
and in the inertial range. Intriguingly, the inertial range signature is not simply that expected 
from homogeneous turbulence, in that there is a robust "envelope" function g(r) for the 
scaling. This is highly suggestive of turbulence in a confined or space-varying medium, see 



for example the work of iBiskamp fc Mullerl (120001 ). Our function g(r) may therefore capture 



the evolution of the turbulence observed at ULYSSES, reflecting both the heating of the fast 
solar wind at the corona and the subsequent expansion in the presence of the large scale 
solar magnetic field. 
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Fig. 1. — Log-log plots of the B-field components' power spectra for days 180 — 189. Two 
regions with different scaling exponents are distingui shable with a break be tween frequencies 
at 10~ 35 — 10 _3 Hz, consistent with previous results (IHorbury et al.lll996al ). For comparison 
purposes, the —1 and —5/3 power scaling laws are also shown. The power spectra for the 
other time intervals examined show similar behavior. 
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Fig. 2. — Log- log plots of third order structure function S3 versus sampling interval r for all 
three components of magnetic field fluctuations in the solar wind measured by the ULYSSES 
spacecraft during contiguous intervals of ten days, which are plotted separately on each panel, 
from day 180 to day 239 of 1995. Only the 0.1% conditioned data is shown for clarity. The 
evolution of the spectral breakpoint to larger r with increasing radial distance (or time) can 
be seen from the change in curvature of logS% between the different time intervals. 
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Fig. 3. — Log-log plots of third order structure function S3 versus sampling interval r 
for all three components of magnetic field fluctuations in the solar wind measured by the 
ULYSSES spacecraft during contiguous intervals of ten days, which are plotted separately 
on each panel, from day 180 to day 239 of 1995. The raw data is represented by "•" 
and the 0.1% conditioned data by "x". Top left panel: radial field Br. Top right panel: 
tangential field B?. Bottom panel: normal field B^. For r > 30 minutes, corresponding to 
power spectral density, this gives evidence for scaling of the form equation (T5]) with 
Cl,b(3) = 0.399 ±0.011, Cl,t(3) = 0.508 ±0.010 and (l,n(3) = 0.445 ±0.008; solid lines show 
linear regression fits for 26 minutes > r > 49 minutes for 0.1% conditioned data, whereas 
the dashed lines show the same fits for the raw data. For r < 30 minutes, corresponding to 
the inertial range, there is no scaling but the data suggest a possible common g{r) defined 
by equation (j3J). 
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Fig. 4. — Evidence for limited variation in ((3) between ten-day data runs in the scaling 
region, and for a common g{r) between the different time intervals. Data for Br, Bt and Bn 
from Figure [2] are replotted after normalization such that the value of S^(t = 30 minutes) is 
the same for each ten-day run within each panel. Gradients in the scaling range (solid lines) 
do not vary systematically with their timing with respect to the 60-day observation period, 
nor do they vary in the same way for different magnetic field components. Away from the 
scaling region, points appear to lie close to a common curve, reflecting g(r). 
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Fig. 5. — Dependence of the value of S^(t = 30 minutes) on sampling time - a proxy for 
mean location of the moving ULYSSES spacecraft - for Br, Bt and B^. Here interval one 
refers to days 180 to 189 of 1995, and interval six to days 230 to 239. 
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Fig. 6. — Evidence for extended self-similarity (ESS) across the full r range. Log-log plots 
of second order structure function S2 versus third order structure function S3 for all three 
components of magnetic field fluctuations in the solar wind measured by the ULYSSES 
spacecraft during contiguous intervals of ten days, which are plotted separately on each 
panel, from day 180 to day 239 of 1995. The different intervals have been uniformly shifted 
in the y-direction for clarity. Top left panel: radial field Br. Top right panel: tangential field 
Bt- Bottom panel: normal field B^. Data points in the inertial range are marked by crosses, 
and in the range by open circles. The straight lines show linear regression fits across 

the full temporal range from r = 2 — 49min. These results imply a global average fitting 
across the different time intervals of Cfl(2)/Cfl(3) = 0.749±0.004, £r(2)/Cr(3) = 0.759±0.004 
and 0v(2)/0v(3) = 0.765 ± 0.004. 
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Fig. 7. — Variations of £(2)/£(3) for fits to the full r range (r = 2 — 49min.) for Br, Bt and 
Bn. Here interval one refers to days 180 to 189 of 1995, and interval six to days 230 to 239. 
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Fig. 8. — Like Figure [6] but showing fits for the inertial range between r = 2 — 14 min. 
These results imply that on average across the different time intervals (s,r(2) / (s,r(3) = 
0.747 ± 0.008, Cs,t(2)/Cs,t(3) = 0.757 ± 0.007 and (sA 2 )/(s,n(3) = 0.753 ± 0.005 in the 
inertial range. 
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Fig. 9. — Like Figure [6] but showing fits for the "/ 1} ' range between r = 26 — 49 min. These 
results imply that across the different time intervals 0.675±0.017 < (l,r(2)/(l,r(3) < 0.745± 
0.017, 0.687±0.019 < Cl,t(2)/Cz,t(3) < 0.835±0.034 and 0.724±0.004 < Cl,tv( 2 )/Cl,jv(3) < 
0.817 ± 0.012 in the "J" 1 " range. 
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Fig. 10. — Variations of £(2)/£(3) for fits to the inertial range (solid lines, r = 2 — 14min.) 
and the range (dashed lines, r = 26 — 49min.) for .Br, St and B^. Here interval one 

refers to days 180 to 189 of 1995, and interval six to days 230 to 239. 
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Fig. 11. — Variations of A for a and (3 quadratic fits (see equations (J7J) and (jHJ)) to the inertial 
range (r = 2 — 14min.) for Br, Bt and B^. Here interval one refers to days 180 to 189 of 
1995, and interval six to days 230 to 239. 
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Fig. 12. — Evidence for g(r) dependence in the inertial range. Log-log plots of the GSFs 
for moments m = — 6 versus g{r) for all three components. A single time interval is 
shown here, days 180 — 189; the same analysis was done for the other time intervals and 
similar results found. The g(r) is normalized for each component such that Cs(3) = 1. This 
is achieved by obtaining Cs(3) from a linear regression fit of the GSF for each component 
and then incorporating this into the g(r) expression such that g(r) = t~ 1 ° 910 ^ Xt ^ 3 ' where 
i = R, T or N. The break between the inertial range and ranges does not seem to 

change position from the GSF but is clearer in these plots. 



